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^ A recent study of dissipative collapse considered a contracting sphere in which the areal 

' and proper radii are equal throughout its evolution. The interior spacetime was matched 

£""""» , to the exterior Vaidya spacetime which generated a temporal evolution equation at the 

00 ■ boundary of the collapsing sphere. We present a solution of the boundary condition which 

£Nj ' allows the study of the gravitational and thermodynamical behaviour of this particular 

* ' radiating model. 
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1. Introduction 

The study of dissipative gravitational collapse achieved prominence with the pre- 
sentation of the junction conditions by Santos^ in 1985. Earlier work on collapsing 
stars in general relativity assumed the exterior spacetime to be empty and as a con- 
sequence, it was required that the pressure at the boundary vanish. Santos provided 
the general junction conditions required for the smooth matching of a spherically 
symmetric, shear-free spacetime to the exterior Vaidya solution across a time- 
like hypersurface. An important consequence of the matching conditions is that 
the pressure on the boundary of the radiating star cannot be zero. It is assumed 
that the interior of the star is radiating energy in the form of a radial heat flux. 
The junction conditions due to Santos rejuvenated the study of gravitational col- 
lapse and the end-states of radiating stars. The simplistic model of Oppenheimer 
and SnydeJ^ has been generalised to include pressure^, anisotropic stresses, elec- 
tromagnetic fielcP^ and the cosmological constantP. These exact models, although 
simplified, give much insight into the dissipative collapse process as well as physical 
characteristics of the radiating star such as its temperature and luminosity. 
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What makes the study of dissipative collapse of stars particularly difficult is 
the solution to the boundary condition representing the conservation of momen- 
tum across the time-like hypersurface. While many exact solutions for shear-free 
radiating spheres have appeared in the recent literature, there are very few models 
that include the effects of shear in the interior of the star. One of the first exact 
models of a shearing radiating star that allowed for an analysis of the gravitational 
and thermodynamical behaviour of the stellar fluid was found by Naidu et aPl 
However, the model was restrictive in the sense that it was acceleration-free, but 
more importantly, the matter variables such as pressure and density became infi- 
nite at the center of the star. It was pointed out that this model could form part 
of a core-envelope model of a radiating star. Further exact shearing solutions were 
obtained by Rajah and Maharaja in which it was assumed that the particle trajec- 
tories within the stellar core were geodesies. An analysis of the temperature profiles 
for these models reveals unphysical behaviour in that the temperatures closer to the 
surface of the star become negative. A recent study of shearing, dissipative collapse 
considered a model of a spherically symmetric matter distribution in which the 
areal radius is equal to the proper radius throughout the stellar evolutiorM These 
so-called Euclidean stars were shown to exhibit very interesting general properties. 
In this paper we present an exact solution to the boundary condition that deter- 
mines the temporal evolution of a Euclidean star. Our solution allows us to study 
the physical and thermodynamical properties of this class of stars even when the 

stellar fluid is far from equilibrium. Since Euclidean stars are not acceleration-free 

[71 

we are able to draw comparisons with the earlier models of Naidu et aV- and Rajah 
and Maharaja. 

The basic structure of our paper is as follows: Following Herrera et al ® we 
present the Einstein field equations for the interior of the stellar fluid in section two. 
In section three we describe the exterior spacetime given by the Vaidya solution 
and we give the main junction conditions required for the smooth matching of the 
interior and exterior spacetimes. In section four we specify the Euclidean condition 
and a particular solution to the boundary condition. In section five we investigate 
the thermodynamical properties of our model and make comparisons to earlier work 
on shearing, radiating models of gravitational collapse. 

2. Shearing spacetimes 

The interior spacetime is described by the general spherically symmetric, shearing 
metric in comoving coordinates 

ds 2 = -A 2 dt 2 + B 2 dr 2 + R 2 (d0 2 + sin 2 0d(f> 2 ) , (1) 

where A = A(t,r), B = B(t,r) and R = R(t,r) are metric functions yet to be 
determined. The matter content for the interior is described by 



Tap = (ju + P±)V a Vp + P±g a p + (P r - P±)XaX/3 + QaVp + qpV a , 



(2) 
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where /i represents the energy density, P r the radial pressure, P± the tangential 
pressure and q a the heat flux vector. The fluid four-velocity V is comoving and is 
given by 

1 



ye 



The heat flow vector assumes the form 

<? Q = (0, g ,0,0), 

since q a V a =0 ensuring radial heat dissipation. We further have 

X a Xa = 1, X a V a = 0. 

The expansion scalar and the fluid four acceleration are given by 

= V.%, a a = Va-pV , 

and the shear tensor by 



V, 



1 



®{g af } + v a v p ). 



(3) 
(4) 
(5) 
(6) 

(7) 



For the comoving line element ([T]) the kinematical quantities take the following 
forms 



(8) 
(9) 

(10) 



where dots and primes denote differentiation with respect to t and r respectively. 
The nonzero components of the Einstein's field equations for the line element <jTJ) 
and the energy-momentum ([2]) are 
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This is an underdetermined system of four coupled partial differential equations in 
seven unknowns viz., A, B, R, P r , P± and q. 



3. Exterior spacetime and junction conditions 

The exterior spacetime is taken to be Vaidya's outgoing solution, given by^ 
2m(v ) 



ds' 



1 



R 



2dvdR + R 2 (d9 2 + sin 2 



(15) 



where m(v) represents the Newtonian mass of the gravitating body as measured 
by an observer at infinity. The necessary conditions for the smooth matching of 
the interior spacetime ((lj to the exterior spacetime tjl 5|) have been extensively 
investigated. We present the main results that are necessary for modeling a radiating 
star. The continuity of the intrinsic and extrinsic curvature components of the 
interior and exterior spacetimes across a time-like boundary are 

2 

/ /."' \ ' 

1 



m(v)i 



(Prh = <?£• 

Relation (fTT|) determines the temporal evolution of the collapsing star. 



(16) 
(17) 



4. Radiating Euclidean stars 

Following Herrera et aW we impose the condition that the areal radius of any 
spherical surface contained within E, with centre placed at the origin, is equal to 
the proper radius from the center through to r = b, the boundary of the star. This 
implies that 

B = R' 



The Einstein field equations (fn j) - (fT4)) reduce to 

R'\ R 

M 

Pr = 



Pi 




(18) 

(19) 
(20) 



(21) 
(22) 
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and the mass function becomes 




m = — — 



The boundary condition (|T7|) yields 

2 
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(23) 



(24) 



valid on r — b. We now have a system of six coupled partial differential equations, 
viz. (fT8 |) -([22 |) and ([24]) . in seven unknowns. 

We focus on (|24|) as a solution of this equation will yield all the relevant kine- 
matical and physical quantities. In doing so, we note that we are requiring (|24j) to 
hold for all r and not just on the boundary r = b. Equation (|24[) is a nonlinear par- 
tial differential equation in A(r, t) and R(r, t). We could analyse it as a quasi-linear 
partial differential equation in A(r, t) only. Then the general solution of ([24| will 
reduce to a general function of the solutions of two ordinary differential equations. 
Unfortunately, these equations are still difficult to solve. As a result, we provide a 
simple solution to ([24]) by setting 



R = aA. 



(25) 



(Note that this closes our system of partial differential equations as we now have 
seven equations in seven unknown functions.) This assumption leads to 

l 2 



R(r,t) = [d(r) e Alt + C 2 (r) e 



where 



Ai = 



1 + V3 



A 2 = 



i - %/3 



(26) 



(27) 



2a 2a 

which is the general solution to the resulting form of (f2"4"|) with the assumption (|25p 
valid for all r. Utilising solution (|2"6")l , the Einstein field equations (fTT]) - (|14p yield 
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4ae- 4A2t (Xxe^rCx + A 2 C 2 



(31) 

where Ai and A 2 are defined in (l27|) . The magnitude of the shear tensor is given by 



(e^C 1 + C 2 ) 3 (e^ i q+C 2 ) 



(32) 



Relation ([3"2l indicates that the shear vanishes when C\{r) oc C 2 (r). In the next 
section we study the thermodynamical properties of our model. In order to ensure 
that the shear remains finite and nonzero for all time we make the following choice 
for our metric function 

R(r, t) = [(a 2 + r 2 )e Alt + (c 2 + r 2 )e A2t ] 2 , (33) 

where a and c are constants and the A^ were defined earlier. 

5. Thermodynamics 

In this section we investigate the evolution of the temperature profile of our model 
within the context of extended irreversible thermodynamics. The causal transport 
equation in the absence of rotation and viscous stress is 

rh a b q b + q a = ~n (h a b \7 b T + Tu a ) (34) 

where h ab = g ab + u a u b projects into the comoving rest space, T is the local equi- 
librium temperature, k (> 0) is the thermal conductivity, and r (> 0) is the relax- 
ational time-scale which gives rise to the causal and stable behaviour of the theory. 
To obtain the noncausal Fourier heat transport equation we set r = in (|34l) . For 
the metric (Q}, equation (|34|) becomes 

T ( qB j +AqB = -^l. (35) 

In order to obtain a physically reasonable stellar model we will adopt the thermo- 
dynamic coefficients for radiative transfer. Hence we are considering the situation 
where energy is transported away from the stellar interior by massless particles, 
moving with long mean free path through matter that is effectively i n hyd ro dy- 
namic equilibrium, and that is dynamically dominant. Govender et a tl0 | H I have 
shown that the choice 

_ f^rr.-a _ (Pi 



k = jT a t c , t c = l^-j T-°, T = IZ-Lj T C , (36) 

is a physically reasonable choice for the thermal conductivity k, the mean collision 
time between massive and massless particles t c and the relaxation time r. The 
quantities a > 0, (3 > and a > are constants. Note that the mean collision time 
decreases with growing temperature as expected except for the special case a = 0, 
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when it is constant. With these assumptions the causal heat transport equation 
(|35l) becomes 

f3(qB)T-° + A(qB) = -a ±—L . (37) 

This equation was comprehensively studied in the noncausal (/3 = 0) case by Govin- 
der and Govender^ as we ll as in specific causal cases. In the noncausal case, (|37f 
can be solved to yield 

(if) 4 "' = [ A 4 - a qB 2 dr + F(t), a^O (38) 

a J 

In (AT) = — — / qB 2 dr + F{t), (7 = 4. (39) 

a J 

where F(t) is a function of integration which is fixed by the surface temperature of 
the star. Note that T corresponds to the noncausal temperature when (3 = 0. For 
a constant mean collision time (c = 0), (|37|) can be integrated to give the causal 
temperature, ie., 



{AT) A = -- 
a 



J A 3 B(qB)dr + J A i qB 2 dr 



F{t) . (40) 



In (|36p we can think of j3 as the 'causality' index, measuring the strength of relax- 
ational effects, with /3 = giving the noncausal case. 
The effective surface temperature of a star is given by 

1 V_L^ 



< rl >==(^JUsJ' (41) 

where L is the luminosity at infinity and <5(> 0) is a constant. The luminosity at 
infinity can be calculated from 

dm /M . 
L 00 = -- r , (42) 
dv 

where m(v) is given in (|2"3"|) . We are in a position to analyse the evolution of the 
temperature in both the causal and noncausal theories. Figure 1 represents the 
causal temperature (dashed line) and noncausal temperature (solid line) as function 
of the radial coordinate. It is clear that the temperature in both the causal and 
noncausal theories is a maximum at the center of the star and drops off smoothly 
as the radial coordinate increases towards the boundary. This trend also indicates 
that the surface layers of the star are much cooler than the interior portions. As 
in the acceleration- free case studied by Naidu et aPand Rajah and Maharaja, the 
causal temperature is everywhere higher than its noncausal counterpart at each 
interior point of the star. The causal and noncausal temperatures are equal at the 
boundary of the star. Figure 1 also reveals that relaxational effects account for a 
larger temperature gradient within the stellar core. This is expected at late times 
during the collapse as the stellar fluid is far from hydrostatic equilibrium. Figure 2 
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illustrates the trend in the relaxation times for the shear stresses. Following Naidu 
et aPJ, the shear transport equation yields 

T i = ir- > (43) 

where coefficient of shear viscosity for a radiative fluid 

V=^r T*T U (44) 
15 

was utilised. In (|43| we have used P = | (Pj_ — P r ) and ro is the radiation constant 
for photons. We have further assumed that t\ — (3\t c . Figure 2 clearly shows that 
the relaxation time for the shear stresses can vary as much as a factor of 10 2 during 
the evolution of the collapsing fluid. A similar result was found for the acceleration- 
free model investigated by Naidu et aP. Figure 3 shows the proper radius as a 
function of time. It is a monotonically decreasing function as expected since the 
star is losing mass in the form of a radial heat flux. It is interesting to note that the 
formation of the horizon can be avoided in our model even in the presence of shear, 
by carefully choosing the arbitrary functions C\(r) and C*2(r). Such a choice would 
ensure that the mass-to-radius ratio, 2ms/rs < 1 which avoids the appearance 
of the horizon for all time. The horizon-free model of a radiating, shear-free star 
undergoing collapse was first studied by Banerjee et oP^. The physical viability 
of this model was studied by Naidu and Govender^l where it was shown that 
the temperature and luminosity profiles were well behaved throughout the stellar 
interior. 

In conclusion, we have presented an exact solution that completely describes 
the temporal and radial behaviour of a particular class of radiating stars, the so- 
called Euclidean stars. We have shown that the model is reasonably well-behaved 
throughout the collapse process, with the physical and thermodynamical variables 
remaining physically viable. Our model of a radiating star with nonvanishing shear 
adds to the very limited class of such solutions that are currently available in the 
literature. 
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Fig. 2. Relaxation time for the shear stress (close to equilibrium - dashed line), (far from equi- 
librium - solid line) versus r. 



September 16, 2010 0:14 WSPC/INSTRUCTION FILE euclideanfinaljunc 



Thermal behaviour of Euclidean stars 11 




8xl0 6 - 



6xl0 6 - 



4xl0 6 - 



2xl0 6 - 



-10 -8 -6 -4 

Fig. 3. Proper radius versus time. 



